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Abstract 



This paper is devoted to the memory of Professor George M. Zaslavsky passed away 
on November 25, 2008. In the held of discrete maps, George M. Zaslavsky introduced a 
dissipative standard map which is called now the Zaslavsky map. G. Zaslavsky initialized 
V— ' \ many fundamental concepts and ideas in the fractional dynamics and kinetics. In this 

^ ■ paper, starting from kicked damped equations with derivatives of non-integer orders we 

derive a fractional generalization of discrete maps. These fractional maps are general- 
izations of the Zaslavsky map and the Henon map. The main property of the fractional 
differential equations and the correspondent fractional maps is a long-term memory and 
^ . dissipation. The memory is realized by the fact that their present state evolution depends 

OO \ on all past states with special forms of weights. 



1 Introduction 



There are a number of distinct areas of mechanics and physics where the basic problems can 
be reduced to the study of simple discrete maps. Discrete maps have been used for the study 
of dynamical problems, possibly as a substitute of differential equations [H El [3l HJ [5] . They 
/\ ' lead to a much simpler formalism, which is particularly useful in computer simulations. In 

c^ ■ this chapter, we consider discrete maps that can be used to study the evolution described by 

fractional differential equations [6l El [9]. 

The treatment of nonlinear dynamics in terms of discrete maps is a very important step in 
understanding the qualitative behavior of continuous systems described by differential equa- 
tions. The derivatives of non-integer orders [6] are a natural generalization of the ordinary 
differentiation of integer order. Note that the continuous limit of discrete systems with power- 
law long-range interactions gives differential equations with derivatives of non-integer orders 
with respect to coordinates (see for example, [10l[TT]). Fractional differentiation with respect 
to time is characterized by long-term memory effects that correspond to intrinsic dissipative 
processes in the physical systems. The memory effects to discrete maps mean that their present 
state evolution depends on all past states. The discrete maps with memory are considered in 
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the papers O [131 [El [IS [lEl [El [TBI [191 120] • The interesting question is a connection of frac- 
tional equations of motion and the discrete maps with memory. This derivation is reahzed for 
universal and standard maps in [181 [E]- 

It is important to derive discrete maps with memory from equations of motion with fractional 
derivatives. It was shown [21] that perturbed by a periodic force, the nonlinear system with 
fractional derivative exhibits a new type of chaotic motion called the fractional chaotic attractor. 
The fractional discrete maps [TSl |T2] can be used to study a new type of attractors that are 
called pseudochaotic [21] . 

In this capter, fractional equations of motion for kicked systems with dissipation are consid- 
ered. Correspondent discrete maps are derived. The fractional generalizations of the Zaslavsky 
map and the Henon map are suggested. 

In Sec. 2, we give a brief review of fractional derivatives to fix notation and provide a conve- 
nient reference. In Sec. 3, the fractional generalizations of the Zaslavsky map are suggested. A 
brief review of well-known discrete maps is considered to fix notations and provide convenient 
references. In Sec. 4, the fractional generalizations of the Henon map are considered. The 
differential equations with derivatives of non-integer orders with respect to time are used to 
derive generalizations of the discrete maps. In Sec. 5, a fractional generalization of differential 
equation in which we use a fractional derivative of the order < /3 < 1 in the kicked term, i.e. 
the term of a periodic sequence of delta-function type pulses (kicks). The other generalization 
is suggested in [18j. The discrete map that corresponds to the suggested fractional equation of 
order < /3 < 1 is derived. This map can be considered as a generalization of universal map 
for the case < /3 < 1. In Sec. 6, a fractional generalization of differential equation for a kicked 
damped rotator is suggested. In this generalization, we use a fractional derivative in the kicked 
damped term, i.e. the term of a periodic sequence of delta-function type pulses (kicks). The 
other generalization is suggested in [18]. The discrete map that corresponds to the suggested 
fractional differential equation is derived. Finally, a short conclusion is given in Sec. 7. 

2 Fractional derivatives 

In this section a brief introduction to fractional derivatives are suggested. Fractional calculus 
is a theory of integrals and derivatives of any arbitrary order. It has a long history from 
1695, when the derivative of order a = 1/2 has been described by Gottfried Leibniz. The 
fractional differentiation and fractional integration goes back to many mathematicians such 
as Leibniz, Liouville, Grunwald, Letnikov, Riemann, Abel, Riesz, Weyl. The integrals and 
derivatives of non-integer order, and the fractional integro-differential equations have found 
many applications in recent studies in theoretical physics, mechanics and applied mathematics. 
There exists the remarkably comprehensive encyclopedic-type monograph by Samko, Kilbus and 
Marichev, which was published in Russian in 1987 and in English in 1993. The works devoted 
substantially to fractional differential equations are the book by Miller and Ross (1993), and the 
book by Podlubny (1999). In 2006 Kilbas, Srivastava and Trujillo published a very important 



and remarkable book, where one can find a modern encyclopedic, detailed and rigorous theory of 
fractional differential equations. The first book devoted exclusively to the fractional dynamics 
and application of fractional calculus to chaos is the book by Zaslavsky published in 2005. 

Let us give a brief review of fractional derivatives to fix notation and provide a convenient 
reference. 

2.1 Fractional Riemann-Liouville derivatives 

Let [a, b] be a finite interval on the real axis R. The fractional Riemann-Liouville derivatives 
D^^ and D'j^_ of order « > are defined p] by 



iD:j)ix) = D:ii:-'^){ 



X] 



1 .o,r\ /M"-.. ,,>„) 



Tin -a) '' .L (x - z)"-"+i 



(D^/)(x) = (-l)"D^(/r")(x) = 



r(n-a) "'i^ (z-x)°-"+i 

where n = [a] + 1 and [a] means the integral part of a. Here D^ is the usual derivative of order 
n. In particular, when a = n & N, then 

(DlDix) = {Dlf){x) = fix), 

{D:J){x) = D^fix), {DU){x) = {-irD-JXx). 

The fractional Riemann-Lioville differentiation of the power functions (x — a)^ and {b — x)^ 
yields power functions of the same form 



a+ 



T{a + /3 + l) 



D^_ib - xf = -fil±^{b - xf-^ 



r(a + /3 + l) 

where P > —1 and a > 0. In particular, if /3 = and a > 0, then the fractional Riemann- 
Liouville derivatives of a constant C are not equal to zero: 



On the other hand, for k = 1,2, ..., [a] + 1, we have 



The equahty 

{D:j)ix) = o 

is vahd if, and only if, 

n 

f{x) = Y,Ck{x-ar-', 

k=l 

where n = [a] + 1 and Ck are real arbitrary constants. The equation 

(D^/)(x)=0 
is satisfied if, and only if, 

n 

f(x) = Y,Ck{b-xr-\ 

k=l 

where n = [a] + 1 and Ck are real arbitrary constants. 

2.2 Fractional Caputo derivative 

The fractional Caputo derivatives '"-D"_|_ and '"-D^. are defined for functions for which the 
Riemann-Liouville derivatives exists. Let a > and let n be given by n = [a] + 1 for a ^ N, 
and n = a foT a & N. li a ^ N, then the fractional Caputo derivatives is defined by the 
equations 



'D:_,f){x) = {i:--D-f){x) = ^ fdz 

V[n - a) Ja 



^x-^l^-'^+i' 






Tin -a) L (z - a;)"-"+i ■ 



where n = [a] + 1. \i a = n E N , then 



( '^D:J){^) = Dlfix), ( ^Z^r-/)(x) = {-iTDlfix). 

\l a ^ N and n = [a] + 1, then fractional Caputo derivatives coincide with the fractional 
Riemann-Liouville derivatives in the following cases: 

( ''d:j){x) = iD:j)ix), 
if 

f{a) = {Dlf){a) = ... = {D--'f){a) = 0, 



and 

( ''DU){x) = W-/)(x), 
if 

fib) = iDifm = ... = {D--'fm = o. 

It can be directly verified that the fractional Caputo differentiation of the power functions 
(x — a)^ and {b — x)^ yields power functions of the form 

^D"" (x - a)^ = ^^^^^^ (x - a)^-" 
'^D^_{b - xf = -i^^±lI-(6 - x)^"^, 



r(a + /3 + l) 

where (3 > —1 and a > 0. In particular, if /3 = and a > 0, then the fractional Caputo 
derivatives of a constant C are equal to zero: 



For k = 0,l,2,...,n-l, 



^D^+C = 0, ^D^^C = 0. 



^D^^{x-a)^ = 0, D^_{b-x)^ = 0. 



The Mittag-Leffler function Ea[X{x — a)°] is invariant with respect to the Caputo derivatives 

""D^+E^Xix - af] = XE4X{x - a)"], 
but it is not the case for the Caputo derivative ^D^_. 

2.3 Fractional Liouville derivative 

Let us define the fractional Liouville derivative on the whole real axis R. The Liouville fractional 
derivatives D'^ and D" of order a > are defined [9j by 

1 9" r f{z)dz 



T{n - a) (9x" J_^ {x - z)°-"+i ' 
(-1)" d"" /•+°° f{z)dz 



T{n - a) ax" J^ {z - x)"-"+i ' 

where n = [a] + 1 and [a] means the integral part of a. Here D^ is the usual derivative of order 
n. In particular, when a = n E N, then 

{Dlf)ix) = iD'_f){x) = f{x), 



(DWix) = D^fix), {Dlf){x) = i-irO^fix). 
If f{x) is an integrable function and f3 > a > 1, then 

{Dll^f){x) = f{x), 

{Dllif){x) = (/r"/)(x), 
{D'Dlf){x) = {±l)\Dl^'^f){x). 
If q; > 0, then the following relations hold: 

{TDif){k) = (Ttkn^fm, 

where 

(=FiA;)° = m" exp{±aTcisgn{x)/2}. 

Here J-" is the Fourier transform. 

2.4 Interpretation of equations with fractional derivatives 

To describe the physical interpretation of equations with fractional derivatives and integrals 
with respect to time, we consider the memory effects and limiting cases widely used in physics: 
(1) the absence of the memory; (2) the complete memory; (3) the power-like memory. 

Let us consider the evolution of a dynamical system in which some quantity F{t) is related 
to another quantity f{t) through a memory function M{t): 

F{t)= f M{t-T)f{T)dT. (1) 

Equation ([1]) means that the value F(t) is related with f(t) by the convolution operation 

F(t) = M(t)*/(t). 

Equation ([1]) is a typical non-Markovian equation obtained in studying the systems coupled to 
an environment, with environmental degrees of freedom being averaged. Let us consider special 
cases of Eq. ([1]) . 

(1) For a system without memory, we have the Markov processes, and the time dependence 
of the memory function is 

M(t-r) = (5(t-r), (2) 

where 6{t — r) is the Dirac delta-function. The absence of the memory means that the function 
F{t) is defined by f{t) at the only instant t. For this limiting case, the system loses all its 
states except for one with infinitely high density. Using ([T]) and ([2]), we have 

F{t)= [ 6{t-T)f{r)dr = f{t). (3) 



The expression ([3]) corresponds to the process with complete absence of memory. This process 
relates all subsequent states to previous states through the single current state at each time t. 

(2) If memory effects are introduced into the system the delta-function turns into some 
function, with the time interval during which f{t) affects on the function F{t). Let M{t) be 
the step function 

M{t-T) =t-\ (0<r<t); 

M(t-r) = 0, (r>t). 
The factor t'^ is chosen to get normalization of the memory function to unity: 

f M{T)dT = 1. 
Jo 

Then in the evolution process the system passes through all states continuously without any 
loss. In this case, 

F{t) = \f^ f{r)dr, 

and this corresponds to complete memory. 

(3) The power-like memory function 

M{t -T)=Mo{t- rf-^ (4) 

indicates the presence of the fractional derivative or integral. Substitution of (j4]) into ([1]) gives 
the temporal fractional integral of order e: 

F{t) = ^j\t-rr-'f{T)dr, (0<^<1), (5) 

where T{e) is the Gamma function, and A = T{e)MQ. The parameter A can be regarded as 
the strength of the perturbation induced by the environment of the system. The physical 
interpretation of the fractional integration is an existence of a memory effect with power-like 
memory function. The memory determines an interval t during which the function /(r) affects 
on the function F{t). 

Equation ([1]) is a special type of equations for F{t) and /(t), where F is directly proportional 
to M * f. In the general case, the values F{t) and f{t) can be related by the equation 

E{F,M*f) = 0, (6) 

where E is a smooth function. Relation ([S]) is a fractional equation for a dynamical system. 
This equation describes the memory effect. If F is a coordinate g(t), either velocity q{t), or 
acceleration q{t), and / is a derivative q^'^\t), then Eq. ([6]) can be considered as an equation 
of motion with memory. For a power-like memory function M{t), we present ([6]) as a equation 

with fractional derivatives: E(q,q,q{t), D'^q) = 0. As the result, we can use the fractional 

calculus [9j to describe the motion of systems. 



2.5 Discrete maps with memory 

The mapping Xn+i = /(a;„) does not have any memory, as the value Xn+i only depends on Xn- 
The introduction of memory means that the discrete value Xn+i is connected with the previous 
values Xn, a;„_i, . . . , Xi. Particularly, any system, which is described by a discrete map, will 
have a full memory, if each state of the system is a simple sum of all previous states: 



^n+l 



E/(^'^)' (7) 



k=l 



where f{x) is a function that define the discrete map. In general, the expression of Eq. ([7]) 
can tend to infinity. Note that the full memory exists for functions that give a finite sum in 
Eq. ([7]). The full memory is ideal because has the same action upon the next states as all the 
others in memory. The map with a long-term memory can be expressed as 

n 

Xn+1 = ^ Vain, k)f{Xk), (8) 

fe=l 

where the weights Va{n, k), and the parameter a characterize the non-ideal memory effects. The 
forms of the functions f{x) and Va{n, k) in Eq. ([8]) are obtained by the differential equation. 
Note that linear differential (kicked) equations gives the linear function f[x) = x. The discrete 
maps with memory is considered, for example, in the papers [121 1131 HSl 111 IISI ^7\ [181 HI 120] . 
The interesting question is a connection of fractional equation of motion and the discrete maps 
with memory fTEj [T^ . It is important to derive discrete maps with memory from equation of 
motion with fractional derivatives. 

3 Fractional Zaslavsky map 

In this section, a brief review of discrete maps is considered to fix notations and provide 
convenient references. For details, see [H [2l [3l HI [5] . Fractional generalizations of these maps 
are discussed. The fractional Zaslavsky maps are suggested. 

3.1 Discrete Chirikov and Zaslavsky maps 

Let us consider the equations of motion 

oo 

i + KGix)Y,s{^-n)=0 (9) 

n=0 

in which perturbation is a periodic sequence of delta-function type pulses (kicks) following with 
period T = 27r/z/, K is an amplitude of the pulses, and G{x) is a some function. This equation 
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can be presented in the Hamiltonian form 

t 



x=p, p + KG{x)J2^{f-^)=0- (10) 

n=0 

It is well-known that these equations can be represented in the form of discrete map (see 
for example Chapter 5 in [2]). The solution of the left side of the n-th kick 

Xn = xitn — 0) = lim xinT — e), Pn = p(tn — 0) = lim p(nT — e), tn = nT (11) 

e-s>0+ e->0+ 

is connected with the solution on the right hand side of the kick x(t„ + 0), p(t„ + 0) by Eq. 
(fTO|) . and the condition of continuity x(t„ + 0) = x(t„ — 0). Using notation (TTTI) . we can derive 
the iteration equations 

Xn+l=Xn+Pn+lT, (12) 

Pn+1 =Pn- KTG{Xn). (13) 

Equations flT2l) and flT3l) are called the universal map. If T = 1 and G{x) = sin(x), then Eqs. 
(11^ and flT^ give the Chirikov map. 

The Zaslavsky map [221 ESI [21 is defined by 

Xn+l =Xn + I^Yn+l + fi, (14) 

r„+i = e-nF„ + £sin(X„)], (15) 

where we use the parameter 

fi=ie^~ l)/q. 

Note that a shift Q does not play an important role and it can be put zero {Q = 0). The 
Zaslavsky map with ^2 = can be represented by the equations 

^n+l = ^n + Pn+1, (16) 

Pn+i = -bPn-Zsm{X^). (17) 

For b = —1 and Z = K, we get the standard map that is described by Eqs. (fT2|) and (fT3|) with 
T = 1. For the parameters 

Z = -£/ie-^ P„ = /iF„, 6 = -e-^ /i = (e« - l)/g (18) 

equations (ITBIl and (TT7I) give Eqs. (TT^ and (TT3]) with fi = 0. For large g — )■ oo (for small 6 — )• 0) 
Eqs. 0161) and (1171) with Z = —K shrink to the one- dimensional sine-map 

Xn+l=Xn + Ksm{Xr,) (19) 

proposed in [25] and studied in many papers. 
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3.2 Fractional universal and Zaslavsky map 

Let us consider a fractional generalization of (Q in the form 



oo 



oD^x + KG{x)J2^{f-n)='^^ (l<a<2), (20) 

n=0 

where o-D° is the Riemann-Liouville fractional derivative [El [HI [9] , which is defined by 

Here we use the notation D^ = d? jdt^, and qI^ is a fractional integration |6llHll9]. The discrete 
map that corresponds to the fractional differential equation of order 1 < a < 2 is derived in 
[mdn]. This map can be considered as a generalization of universal map for the case 1 < a < 2. 



Theorem 1. Fractional differential equation of kicked system ^2D\) the initial conditions 
is equivalent to the discrete map 

Pn+l=Pn-KTG{Xn), (22) 



Xn+1 



^Pfc+iK.(n- fc + 1) + .^— ^(n + l)"-2, (23) 



r(«)^"^"^ "^ ' ' ' r(«-i) 



where 1 < a < 2 and the function Vai^z) is defined by 

V^{z) = z''''-{z-ir-' {z>0). (24) 



Proof of this Proposition is realized in [TS] for C = 0, and in [T^ for C j^ 0. Equations 
and fl23l) define a fractional universal map in the phase space {xn,Pn), which is a generalization 
of the map ([12]) and ([I3]). Note that Eqs. ([22]), ([23]) with G{x) = sin(x) give 

Pn+i = Pn - KT sin(a;„), (1 < a < 2), (25) 

rpa—l " C^l^° 



Xn+1 



5^P.+iK.(n- A; + 1) + _ (n + ir-\ (26) 



Via) ^^^^' "^ ■ ' ■ Via- I] 



These equations define a fractional standard map on the phase space, which can be called the 
fractional Chirikov- Taylor map. 
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The fractional generalization of Zaslavsky map [22J [23] can be defined by 

Pn+i = -bpn - Z sin(x„), (1 < a < 2), (27) 



Xn+i = Yf^ X^Pfc+iV;(n - A; + 1) + _ (n + l)"-^, (28) 

where the parameters are defined by conditions (ITSl) . For b = —1 and Z = K, Eqs. ( l27l) and 
(EH]) give the fractional standard map that is described by Eqs. ( I25l) and ( l26l) with T = 1. 
For for small 6—7-0 and C = 0, the equations (I27p and (I28p with Z = — K shrink to a 
one- dimensional map 

Xn+i = YT^y2Va{n-k + l)sm{xn), (K a < 2). (29) 

1 [a) ^ — ^ 
^ ^ fc=0 

This map can be considered as a fractional generalization of one-dimensional sine- map (I19p . 
The fractional sine-map f l29|) is characterized by a long-term memory such that the present 
state evolution depends on all past states with special forms of weights function ((24 



3.3 Kicked damped rotator map 

Let us consider a kicked damped rotator [4j. The equation of motion for this rotator is 



oo 

+ qx = KG{x)J2^(t-^T\ (30) 



It is well-known that Eq. ( 130|) gives [4j the two-dimensional map 

p„+i = e-«^[p„ + KG(x„)], (31) 

1 - e"'?^ 

Xn+l=Xn-\ [pn + KG{Xn)]. (32) 

Q 
This map is known as the kicked damped rotator map. The phase volume shrinks each time step 
by a factor exp(— g). The map is defined by two important parameters, dissipation constant q 
and force amplitude K. These equations can be rewritten in the form 

p„+i = e-«^[p„ + KG(x„)], (33) 

Xn+l = Xn^ Pn+1- (34) 

It is easy to see that these equations give the Zaslavsky map f fT6|) and f fT7|) with f2 = if 
Xn = Xn, Yn = Pn, s = K, T = 1, G{x) = sin(x). 
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If we consider the limit g — )■ cxd and K ^ oo such that q/ K = 1 and the function 

G{x) = (r — l)a; — rx^, 
then equations (1521) and (15^ give the logistic map 

Note that the evolution of logistic maps with a long-term memory is considered in [IT] . 

3.4 Fractional Zaslavsky map from fractional differential equations 

Let us consider a fractional generalization of differential equation for a kicked damped rotator 
()30p in the form 

oo 

oD^x-qoD^x = KG{x)J2^it-^T), (35) 

n=0 

where q & R and o-D" is derivative of non-integer order. This equation will be considered with 

l<a<2, /3 = a-l, 

and o-D" is the Riemann-Liouville fractional derivative |S1 El E] defined by f l2ip . Note that we 
use the minus in the left hand side of Eq. (135|) . where q can be positive and negative value. 
The discrete map that corresponds to the fractional differential equation has been derived in 
[18]. 



Theorem 2. Fractional differential equation of kicked system [JB^) is equivalent to the discrete 
map 

p„+i = e«^[p„ + irG(xO], (36) 



Via 

^ ' fc=0 



1 '^ 

V pfc+1 iy„(g, T, n + 1 - /c), (37) 

where the functions Wq{(1-,T^ z) are defined by 

iy„(g,T,m + l) =T"-^ /" e-'^^^ {m + y)"-'^ dy . (38) 

Jo 



Proof of this statement is realized in [18] . Equations ( l36l) and ( 1371) define a fractional kicked 
damped rotator map. 
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This Theorem can be used to derive a fractional generahzation of the Zaslavsky map. If we 
use the conditions 

X„ = x„, Yn = Pn, e = K, T = 1, G{x) =sm{x), (39) 

then Eqs. (136|) and (1371) give the map 

F„+i=e'^k + £sin(XOl, (40) 



1 

These equations can be considered as a second fractional generalization of the Zaslavsky map 
f lT^ and flTSj) with fi = 0. Note that this fractional Zaslavsky map is derived from fractional 
nonlinear differential equation f l35|) with condition ( 139|) . If we use the variables 

then discrete equations 



X, 



n+l 



P„+i = -6P„-Zsin(X„), (42) 

1 r). 



-- ^ Pfe+i W^{q, l,n + l-k). (43) 



r(a _, 

^ ' k=0 



These equations can be considered as a fractional generalization of the Zaslavsky map equations 
(IT6|) and flT7|) with f] = 0. For a = 2, this fractional Zaslavsky map gives the Zaslavsky map 
that is described by Eqs. flT6l) . and flTTl) . For exp{g} — )■ equations fHOl) and PUJ) shrink to a 
one-dimensional map that can be considered as a fractional generalization of one-dimensional 
sine-map flT^ . 

4 Fractional Henon map 

4.1 Henon map 

The Henon map [26l [27] can be considered as a two-dimensional generalization of the logistic 
map [1]. Then the Henon map is presented as a map given by the coupled equations 

a;„+i = 1 - axl + i/n, (44) 

l/„+i = 6x„. (45) 

Substitution of Eq. (^ into Eq. (jBl) gives 

x„+i = 1 - aa;^ + bxn-i- (46) 
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This equation is equivalent to the Henon map that is defined by Eqs. ( I44l) and (l45l) . Note that 
the Henon map in the form (H6ll can be derived from the generahzed dissipative map. 
Let us consider the generahzed dissipative map 

Xn+l =Xn+ Pn+lT, (47) 

p„+i = -bpn + KTGixn), \b\ < 1. (48) 

Using Eqs. (H7|) and (HHj) with K = T = 1, and the function 

G{x) = 1 - (1 - 6)x - ax^ (49) 

the Henon map in the form (H6l) can be obtained. Let us prove this statement. Substitution of 
(H7|) in the form Pn+i = Xn+i — Xn into (HH!) gives 

Xn+l - Xn = -&(Xn - X„_i) + G(Xn). (50) 

Then we have 

Xn+l - bXn-i = (1 - b)Xn + G(x„). (51) 

Substitution of fH9|) into flSTj) gives the Henon map P6|) . 



4.2 Fractional Henon map 

To derive a fractional generalization of the Henon map we can use the fractional generalized 
dissipative map and the fact that the Henon map in the form (H6!) can be derived from the 
generalized dissipative map. 

Let us consider a fractional generalization of the map 0471) and (H51) in the form 

rpa — l ^ 

Xn+i = -— -Vpfc+iV;(n-A; + l), (1< a < 2), (52) 

Fa ^-^ 

^ ' fc=0 

Pn+i = -bpn + KTG{xn). (53) 

This is the fractional generalized dissipative map. If a = 2, then Eqs. fl52l) and fl53|) give fHTl) 
and (gHD. If G'(x) = sin(x), T = 1 and /sT = -Z, C = 0, then Eqs. ([52]) and ([53]) give the 
fractional Zaslavsky map equations (|27|) and fl28l) . 

To derive a fractional generalization of the Henon map, we can use the fractional dissipative 
map ( l52l) and ( l53l) with i^ = T = 1, and the function ( H9l) . As a result, we obtain the equations 



1 " 
Xn+l = :p^^y'pfc+i14(n- A; + l), (54) 

1 (a) ^-^ 

^ ' fc=0 

Pn+i = -&Pn + 1 - (1 - &)2;n - cta;^ (1^1 <1), (55) 
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where 1 < a < 2 and Va{z) is defined by (l24l) . These equations can be considered as a first 
fractional generahzation of the Henon map (jSj), (l45l) . 

Fractional Henon map flM|) and fl55|) is derived from Eq. fl52l) and fl53l) . where the dissipation 
has been introduced into the iteration equation. In this case the fractional equation of motion 
with dissipation is not used. Equations f l5^ and f l55|) are not derived from fractional differential 
equations of kicked damped system. We can suggest the second possible getting of a fractional 
generalization of the Henon map. This map can be derived from the fractional equation of 
kicked damped rotator (!35|) . 

It is known (see Sec. 1.2 of [4J) that the Henon map can be derived from equation of motion 
([30D with 

G{x) = -T^ [1 + (1 + b)x + ax^] , (56) 

where b = — exp(— g) and T = 1. Using Eq. (!30l) with function (|56ll . we obtain the map 
equation 

a;„+i = a;„ H J9„ - 1 - (1 + 6)x„ - axl, (57) 

q 

Pn+i = e^V + -^ [1 + (1 + b)x^ + axl] , (58) 

1 + 

where T = 1, and b = — exp{— g}. These equations are equivalent to the Henon map equation 
(jlGjl . Proof of this proposition is realized in |3]. 

Let us consider the fractional differential equation (!35|) that are generalization of Eq. ( |30|1 
by non-integer order derivatives. In order to derive fractional generalization of Henon map, we 
use Eq. 1^^ with the function S5E\) . As a result, we obtain the equations 

1 " 

Xn+i = TT7 TTy^^Pk+iWa{q,T,n + l-k), (59) 

1 a — 1 ^^ 



Pn+i = e'^ 



fc=0 



Q 



Pn + -^ {I + {1 - b)xn + axl) , (60) 

1 + J 

where the functions Wa{q,T,z) are defined by Eq. ( 138|) and q = ln(— 6). These equations 
define a second fractional generalization of the Henon map. This map is derived from nonlinear 
differential equations with derivatives of non-integer order. Equations (1591) and (!60l) can be 
called fractional Henon map. 

5 Fractional derivative in the kicked term and Zaslavsky 
map 

In this section, a fractional generalization of differential equation (jH]) is suggested. In this 
generalization, we use a fractional derivative of the order < /3 < 1 in the kicked term, i.e. 
the term of a periodic sequence of delta-function type pulses (kicks). The other generalization 
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that is described by Eq. (1201) is suggested in [18]. The discrete map that corresponds to the 
suggested fractional equation of order < /3 < 1 is derived. This map can be considered as a 
generahzation of universal map for the case < /3 < 1. 

5.1 Fractional equation and discrete map 

Let us consider a fractional generalization of (Q in the form 

D',x + KG{^D^x)Y,s[j;-n)=0, (0 < /3 < 1). (61) 

n=0 

where qD^ is the Caputo fractional derivative 

o^Dfx = oiyolx = —1- r -^^, (0 < /3 < 1). (62) 



Here we use the notation Dj = d/dt, and qI^ is a fractional integration |6l El |9] . For (3 = 
fractional equation (16T!) gives Eq. IQ. 

Theorem 3. Fractional differential equation of kicked system l[61\) is equivalent to the discrete 
map 

Xn+l = Xn+ Pn+lT, (63) 

Pn+1 =Pn- KTG( _ ^Pfc+iV-2-^(n - A:)J, (0 < /3 < 1), (64) 

'^ ^^ k=0 

where the function V2-p{z) is defined by 

V2-^{z) = z'-^-{z-lY-^ {z>l). (65) 



Proof. The fractional equation f l6T|) can be presented in the Hamiltonian form 

X = p, 

t 



p + KG{ ^D^x) E Kt " '') = °- ^^^^ 

n=0 

Between any two kicks there is a free motion 

p = const, X = pt + const. (67) 
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The integration of (l66l) over (t„ + e, t„+i — e) gives 

x{tn+i - 0) = a;(t„ + 0) - Pn+iT, (68) 

p(t„+i-0)=p(t, + 0). (69) 

The solution of the left side of the n-th kick 

Xn = xitn — 0) = limx(?7,T — e), (70) 

£-5-0 

Pn =p(in - 0) = limp(nT -e), tn = nT, (71) 

is connected with the solution on the right hand side of the kick x(t„ + 0), p(t„ + 0) by Eq. 
(|66|) . and the continuity condition 



x{tn + 0)=x{t„,-0). (72) 

The integration of fl66|) over the interval (t„ — £, t„ + e) gives 

p(t„ + 0) = p(t„ - 0) - KTGi^Dix). (73) 

Using notations (170]) and f lTT]) . and the solution f E7|) . we get 

p(t„ + 0) =p(t„+i-0) =p„+i. (74) 

Substituting of (!72|) and (!73|l into (168|) and (!69l) . we derive the iteration equations 

Xn+l = Xn -\- Pn+l-L ; (75) 

Pn+i=Pn-KTGi^'Dtx). (76) 

To derive a map, we should express the fractional derivative 

through the variables (1701) and (17T!) . Using a; = p, we have 



Cn/3^_ i / " p(r)rfr 



°^*"'^"r(l-/3)yo (tn-rr ^^^ 



It can be presented as 



n-l 



Cn/3 . _ 1 Y- n' Pi^)dr 



where t^+i = ^fc + 7" = (/c + 1)T, tk = kT, and to = 0. 
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For the interval (tfc,4+i), Eqs. (ETj), ([70]) and (EI]) give 

p{t) = p(tk + 0) = p{tk+i - 0) = pk+i, T e (tfc, tfe+i). 



^0) 



Then 



.^ _^0 = Pk+l / (in -t) ^dT 



tn-tk ^1-/3 

Pfc+1 / Z-^dz=Pk+l-—-r 

i„-ti,^i J^ P 



n — tfc + l 



tn—tk 



*n— tfc + 1 



Pfc+1 



1-/3 

2^1-/3 



(tn — tk) — (tn — tk+l) 



Pk+l 



1-/3 



{n-ky-" -{n-k-1) 



1-/3 



(81) 



Using (1 — /3)r(l — /3) = r(2 — /3) and Eq. f lHTj) . the fractional derivative f l79|) can be presented 
as 



A„^ 



yi 



r(2-/3) 



5^Pfc+iV2-/3(n-A;), (0</3<l), 



A:=0 



where 



V2-f,{z)=z'-^-{z-lY-^. 



(821 



^3) 



This ends of the proof. D 



As a result, Eq. ( !76|) takes the form of Eqs. (163|) and (IMI) . These equations define the 
fractional generalization of universal map. Let us consider the following Corollary. 

Theorem 4. For /3 = the fractional universal map, which is defined by Eqs. (E^j and ( [g^ , 
g'Zfes i/ie usual universal map. 

Proof. Let us consider that the fractional universal map (p5]) and (^^ for /3 = 0. Substi- 
tution of Eq. f l63|l in the form 

Pk+l = j;{xk+i - Xk) (84) 



into iteration equation (j64ll gives 

p/2_ /g) XI (^'^+1 ~ Xk)V2-p{n - k) 



Pn+l 



n-l 



(85) 



A:=0 



Then, the fractional map ( l63l) . ( I64l) is defined by 

Xn+l = Xn -\- Pn+lJ- i 
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Pn+i =Pn-KTG(^———^{xk-Xk^i)V2-(3{n-k)j. (86) 

where < /3 < 1. For /3 = 0, we have V2-j3{z) = 1, and 

n-l 

^(Xfc+i - Xfc) = Xn- Xq. (87) 

fc=0 

Then Eq. (!85|) gives 



J5„+l = Pn - -f^T G f x„ - Xo j . 



As a result, Eqs. fl86|) with /3 = give the usual map for the case Xq = 0. 
This ends of the proof. D 



5.2 Examples 

Let us consider some examples of the fractional map (163|) . ( 164|) . 

Example 1. If G(x) = -x, then Eqs. (ES]), (EH) are 

^^n+l = 3;n + Pn+lT, 
rp2-l3 "--^ 



„+i = p„ + K:^—-—J2Pk+iV2^^{n -k), {0</3< 1), (89) 



where ^2-/3(2;) is defined in fl24|) . Equations (!89|) can be considered as a fractional generalization 
of the Anosov type system. 

Example 2. If G{x) = sin(a;), Eqs. (ES]), ([MD give 

Xn+l = Xn -\- Pn+lT, 

rpl-13 "Z^ 

Pn+1 =Pn + KT smi^—-—Y,Pk+iV2-p{n -k)j, {0 < f3 < 1). (90) 

This map can be considered as a fractional generalization of standard map. The other possible 
form of Eq. (!90|) is 

-. 71—1 

p„+i = p. - i^ sin( ^ 5^(2:^+1 - Xfc) [(n - ky-^ ~{n-k- 1)^-^] ) , (91) 



-r(2-/3),.o 
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where we use Eqs. (18^ . (12^ . T = 1, and < /3 < 1. These equations define a fractional 
standard map that can be called the fractional Chirikov map. 



Example 3. The fractional generalization of Zaslavsky map [221 123] can be defined by 

Xn+l = Xn + Pn+l, (92) 

Pn+1 = -bpn- 



— K sin 



(E%fr|r[(--^)"'-(--^-i)"1)- (93) 



fc=0 



For 6 = 1, we get the fractional standard map fPT]) . This map is one of possible fractional 
generalizations of the Zaslavsky map. In this generalization we introduce a dissipation by the 
change of the variable pn — > —bpn- Equation (193!) is not directly connected with a fractional 
equation of motion. A generalization of the Zaslavsky map that is derived from the differential 
equation with fractional damped kicks is suggested in the next section. 

6 Fractional derivative in the kicked damped term and 
generalizations of the Zaslavsky and Henon maps 

In this section, a fractional generalization of differential equation fl5Ul) for a kicked damped 
rotator is suggested. In this generalization, we use a fractional derivative in the kicked damped 
term, i.e. the term of a periodic sequence of delta-function type pulses (kicks). The other 
generalization, which is described by Eq. (!35|) . is suggested in [18]. The discrete map that 
corresponds to the suggested fractional differential equation is derived. 

6.1 Fractional equation and discrete map 

Let us consider the fractional generalization of equation fl30l) in the form 



oo 



D^x - qD]x = KG{ qD^x) ^ 6{t -nT), (0 < /3 < 1), (94) 

n=0 

where q E R, and qD^ is the Caputo fractional derivative [H] of the order < /3 < 1 defined 
by (IS^ . Note that we use the minus in the left hand side of Eq. fl5^ . where q can be posi- 
tive and negative value. The fractional deribative qD^x is presented in the kicked damped term. 



Theorem 5. Fractional differential equation of kicked system (f^l) is equivalent to the discrete 
map 

l-e-"^ . X 

Xn+l =Xn^ Pn+l, (95) 

q 
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n-1 



Pn+i = e 



qT 



Pn + KG(———'^Pk+iW2-p{q,T,n-k) 



r(i-/3) 



k=0 



where the functions W2-B o?^e defined by 



W2-^{q, T, m) = T^-^ f e^^^^"^) (m - z)"'' dz. 

Jo 



(96) 



(97) 



Proof. Fractional equation flMjl can be presented in the Hamiltonian form 

X = p, 



p-qp = KG{^D^x)J2^it-nT), (0 < /3 < 0, gG 
Between any two kicks 



n=0 



p — gp = 0. 
For t G (tn + 0, t„+i - 0), the solution of Eq. ([99]) is 

p{tn+i - 0) = p(t„ + 0)e'?^. 

Let us use the notations tn = nT, and 

Xn = x{tn — 0) = lim x(nT — e), 

Pn = p(tn - 0) = lira p{nT - e). 

For t E {tn — £,tn+i — e), the general solution of (p8l) is 

00 „t 

p(t) = p^e''^*-*") + ir J] G'( ?Aia;) / rfre''(*-")<5(r - mT). 

m=0 >^i7i-e 



Then 



0IT 



Ct^H 



Pn+i = e^^ Pn + KG{ A„a; 
Using (I103P , the integration of the first equation of (!98|) gives 

1 - e^^ r 



Xn-|_i Xr, 



q 



Pn + KG{^Dlx 



Let us consider the Caputo fractional derivative from Eqs. (11031) . (I104p . Then 

1 /■*" dr dx{T) 



^dIx =0 ll-^Dlx 



r(l-/3)io {tn-rf dr 
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-, (0</3<l) 



(98) 

(99) 

(100) 



(101) 



(102) 



(103) 



(104) 



rl-/3. 



Using a; = p, we have 

^t„^ -0 -'t ^Yii 
This relation can be rewritten as 



p{T)dr 
(l-/3)7o Jt~^W- 



"oDtx 



^ n— 1 



tk 



+^ p{T)dT 



where tk+i = tk + T = {k + 1)T, and 4 = /cT, such that to = 0. 
For r e (tfc,tfc+i), Eqs. (dUQ]) and (^2) give 



Pfc+ie''(^-*'=-^)=Pfc+ie''(^-*'=+^). 



Then 



*fc+i 



'tn—tk + l 



tk + 1 



p{T)dT 



g{T-tk+i)(f _ 



(t„ - T)-"dT 



litn-tk+l-z)^-l3^^ ^ p^^^^ 



tn—tk 

q{tn-tk+i) I z'/^e-^^dz 

tn—tk + l 



Pk^ir -'e'^^-'-'^^ / y-^e-ydy. 



g(tn— tfc + l) 



As a resuh, Eq. fll07p can be represented by 



■k+ 



Pfc+ig^-ie''^"-^-!)^ 



^ p{T)dT 
tk K^^^ 

r(l - /3, g(t„ - tfc+i)) - r(l - /3, g(t„, - tk)) 



(105) 



(106) 



(107) 



:io8) 



Here T{a, b) is the incomplete Gamma function, where a, b are complex numbers. This function 

can be defined by 

b"- 

T{a, b) = T{a) iFi(l, 1 + a; -b). 

a 

Here iFi is the confluent hypergeometric Kummer function |28j . 

{a)kz'' 



iFi{a,c;z) = ^ 



fc=0 



'c)k k\ 



and {a)k is the Pochhammer symbol 

{a)k = a{a + l)...{a + k -1), k e N. 
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Using fITOHll and fITnSll . we obtain 



^ 71—1 

K^ = p7Y3^$^Pfc+iW^2-/3(g,T,n - /c), (0 < /3 < 1), 







where PF2_/3(g, T, m) is defined by ([9] 

Substitution of f lT09|) into flT03|) and fITOil) gives 



n-l 



Pn+1 = e 



qT 



Pn + KG(^^r^— ^ ^Pfc+iVr2-/3(g, T,n- k)j 



k=0 



l-e-?^ 



n-l 



•^n+1 -^n 



Pr, 



'n + KG[—^-—Y,Pk+iW2-p{q, T,n-k) 



k=0 



Equations f lllOp and (11 111) can be presented in the form of Eqs. fl95|) and fl95|) . 
This ends of the proof. D 



(109) 



(110) 



:iii) 



6.2 Fractional Zaslavsky and Henon maps 

The iteration equations fl^S]) and (JHSD define a fractional generahzation of the kicked damped 
rotator map fl5(I]) . If we use the conditions fl5^ . then Eqs. fl^ and flU^ give the map 






Yn+1 



Yn + e sin 



1 



where 



/^ = 



Y,yk+iW2-(,{q,T,n-k) 



(112) 
(113) 



These equations can be considered as a fractional generalization of the Zaslavsky map ([Til) 
and flTSl) with = and fi = fi'. This fractional Zaslavsky map is derived from fractional 
differential equation flMj) with condition (|39|) . For (3 = this map gives the Zaslavsky map (IT^ 
and ( fTSl) with n = jj! and fi = 0. 

By analogy with Sec. 3, we can derive a fractional generalization of the Henon map from 
equation of motion ( 194]) with 



G{x) 



1 + 6 



[l + (1 + 6)a; + ax^] 



(114) 



where 6 = — exp(— g) and T = 1. As a result, this generalization has the form of Eqs. f l95|l 
and fl96l) with the function (I114p . This fractional Henon map is derived from the differential 
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equations with fractional derivative in the kicked damped term. For /3 = this map gives the 
usual Henon map. 

Computer simulations of the fractional generalizations of discrete maps prove that the non- 
linear dynamical systems, which are described by the equations with fractional derivatives, 
exhibit a new type of chaotic motion. This type of motion can be considered as a fractional 
generalization of chaotic attractor. As a result, the fractional discrete maps allow to study a 
new type of attractors that are called pseudochaotic. 

7 Conclusion 

In many areas of mechanics and physics the problems can be reduced to the study of discrete 
maps. In particular the special case of discrete maps has been studied to describe properties 
of regular and strange attractors. Under a wide range of circumstances such maps give rise 
to chaotic behavior. The suggested fractional maps can be considered as a fractional gener- 
alization of discrete map. Note that a wide class of these maps can be derived from kicked 
fractional differential equations. The suggested fractional Zaslavsky map and the fractional 
Henon map can demonstrate a chaotic behavior with a new type of attractors. The interesting 
property of these fractional discrete maps is a long-term memory. As a result, a present state 
evolution depends on all past states with the fractional power-law weights functions Va{z) and 
Wa{q,T, z) defined by equations fl2^ and fl38|) . Note that the fractional maps are equivalent 
to the correspondent fractional kicked differential equations. An approximation for fractional 
derivatives of these equations is not used. This fact can be used to study the evolution that is 
described by fractional differential equations and to describe pseudochaotic attractors. 
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